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NOMENCLATURE 


A Local cross-sectional flow area* pre-exponential term 

in Arrhenius equation 

A Q cons tant 

bi»b 2 blowing parameters 

C* characteristic velocity 

C^- skin friction coefficient 

Cp specific heat at constant pressure 

D local diameter* damping function 

d penetration depth 

E/R characteristic activation temperature 

h enthalpy* heat transfer coefficient 

k thermal conductivity 

K m mixing length constant 

m exponent of temperature dependence of viscosity 

M Mach number 

u density ratio exponent, Arrhenius equation 

p pressure 

p + pressure gradient parameters - 

P„ Prandtl number 

r 

q heat flux rate 

Q heat of pyrolysis 

r„ radius of curvature at the throat 

c 

r nozzle radius, at any section 

recovery factor 



R 


Stanton number 


’w 


(ht.adV ?e u « 


time 

temperature 

dimensionless velocity u/u^ 

velocity component tangential to the wall 

friction velocity ( 3 w )^ 

velocity component normal to the wall 

distance along the wall 

Char depth 

distance normal to the wall 

$w u x y 

dimensionless normal distance 

HW 


axial distance from throat 
nozzle half angle 
ratio of specific heats 


70 


displacement thickness v* 


<x> 


momentum thickness 


Su 
? e e 


o 


(i- ^r><»y 
> 0 0 


(l-~7p) dy 
e 


thermal diffusivity 
dynamic viscosity 
kinematic viscosity 
density 
shear stress 
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T 


CO 




dimensionless velocity = u/u 


c<5 


energy thickness ^ 


3 u 


*> -u 


(1- 


e e 


viscous transformation exponent 

h, - h 


dimensionless enthalpy 


t ,w 


h t,o“ h t,w 


Subscripts 


ad 

c 

e 

1 

o 

P 

w 

t 


adiabatic wall 
char 

free stream 

virgin material 

stagnation condition 

pyrolysis 

wall 

total 


Superscripts 

o initial conditions 

* 


^)<*y 


throat conditions 



ABSTRACT 


The determination of heat transfer and char layer 
thickness in a composite liner of a rocket nozzle has been 
done by using approximate methods. The compressible 
turbulent boundary layer flow along the nozzle wall is 
analyzed to obtain heat transfer rate by using modified 
Van Driest* s eddy diffusivity model for velocity profi^,. 
quadratic relationship between the velocity and the total 
enthalpy as well as the integral momentum and enthalpy 
equations. Assuming quasi-steady situation the char layer 
thickness is calculated by using a simplified procedure 
for approximate solution of heat transfer problem with 
a moving interface. Comparison of the calculated results 
with some experimental data and other related theoretical 
results has been done. 



CHAPTER 1 


INTRODUCTION 

1.1 GENERAL BACKGROUND 

During the operation of a chemical rocket motor, heat 
is transmitted to all the parts which are exposed to the 
combustion products, such as, solid-propellant grain, 
liquid-propellant injector face, chamber and nozzle 
walls. The determination of heat transfer rate and the 
choice of an appropriate thermal protection method is 
essential for the design and development of rocket motors. 
Here we will be concerned, only with the heat transfer 
and cooling in a rocket nozzle. The objective is to 
protect the nozzle walls especially the critical hot 
regions so that the performance of the rocket motor 
follows the desired behaviour. 

Inadequate thermal protection of a nozzle may lead 
to the breakdown/unpredic table rocket motor operation, such 
as, in the following situations ; 

1. The wall temperature exceeds the value at which the 
material is readily melted or oxidized. The local 
loss of material weakens the wall so that it is no 
longer able to withstand the imposed load. 
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2, An uneven erosion of the throat region takes place 

resulting in mis-alignment of the thrust vector. 

Some of the cooling methods used in rocket nozzles are j 
1) Insulation cooling, 2) Heat sink cooling, 3) Transpiration 
cooling, 4) Regenerative cooling, and 5) Ablative cooling. 
Some times a combination of some of these methods has also 
been employed. 

For solid propellant rocket motors, only a few of the 
above-mentioned methods like insulation cooling and ablative 
cooling have been found useful. The ablative cooling has 
been successfully used in large booster rockets and the 
present study deals only in this method of cooling. 

In ablative cooling, the inner surface of the nozzle 
walls is lined with a composite material consisting of 
layers of woven fibres (such as silica or carbon/graphite 
fibres) impregnated with an organic plastic material (such 
as epoxy or phenolic). The liner when exposed to a high 
heat flux undergoes complex thermo-physical phenomena 
keeping the inner nozzle wall at the desired lower temp 
The term ’ablation’ implies the occurrence of heat and 
mass transfer processes in which a large amount of thermal 
energy is expended by the loss of surface region material. 

A simplified schematic diagram identifying some distinguishing 
features for an ablative liner is shown in Fig. l. 
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When the ablator surface temp, rises, the heat is 
transmitted into the virgin material and the organic plastic 
component first begins to decompose (pyrolysis zone II) 
producing gaseous products which diffuse outwards through the 
matrix of higher temp resistant fibres forming the so-called 
char layer (zone III \ As the surface temp increases further, 
other chemical reactions such as evaporation or oxidation of 
char material set in, resulting in removal or erosion of 
surface material (zone IV), Aerodynamic forces as well as 
particle impact also contribute towards the surface erosion. 
Evidently the rate of char layer formation or surface regre- 
ssion depends upon the temp .velocity, etc. of the flow 
field over the ablator surface, and consequently, the 
ablator requirements will vary from one situation to another. 
In other words, large erosion rates in reentry heat shields 
may be desirable but not so in the case of rocket nozzle. 

Here the analysis is restricted to the determination of char 
layer thickness for a given heat transfer rate by using a 
highly simplified approach. Before giving the details of 
the present analysis, the past literature relevant to the 
topic is reviewed next. 

1.2 REVIEW OF PAST CONTRIBUTIONS 

The phenomena associated with the ablative cooling in 
rocketry has been extensively studied for the past two 
decades. The analysis of the problem involves determination 
of the heat flux at the liner surface exposed to hot 
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gases flowing over it and of the effect of various thermo- 
physical process occurring inside the liner material. Most 

of the investigations usually deal with only one part of 
the problem, that is, either the heat transfer rate or the 
degradation of the liner, in details by making use of the 
available experimental or theoretical results for the other 
part. 


Some of the earliest investigations for the prediction 
of heat transfer rate from the combustion products to the 
rocket nozzle walls have been done and reviewed by Bartz [l] 
Using integral method and l/7th power dis tributions of 

velocity and temperature in terms of certain boundary 
layer shape parameters, Bartz obtained the following 
expression for the local heat transfer coefficient. 


h (; r &^ 0 - 8 <£ )0 * 1]{ £ >0,9 

r w 


( 1 . 1 ) 


where 


a = 


[«£>C1+ m2 )+ M 0 - 8 -WV [1 + 3^1 


The evaluation of the constant 0.026 as well as the 
dependence on the wall temperature was achieved by comparing 
the calculated results with the available experimental data 
for both heated air and rocket thrust chambers. Bartz 
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expression has been widely used because of its simplicity 
and has been in general found to over-estimate the heating 
rate. 

During the past decade or so, substantial progress 
has been achieved in the description of compressible turbulent 
boundary layer flows with pressure gradients and heat as well 
as mass transfer due to improved experimental and theoretical/ 
numerical techniques [2,3]. More elaborate methods of 
predicting heat transfer rates make use of the semi- 
empirical models for turbulent transport coefficients 
appearing in governing differential conservation equations. 

The differential methods involve excessive amount of numerical 
computation on a digital computer, and complexity of the 
solution increases enormously. When coupled with the 
phenomena occurring in the liner. It is therefore found 
worthwhile to devise approximate methods for predicting 
rapidly but reasonably correct heat transfer rates into 
converging-diverging nozzles. 

Recently, Mastanaiah [4] reported an approximate 
method based on the integral equations and the eddy diffu- 
sivity model due to Van Driest and Cebecci [3] and it is 
claimed that the heat transfer results are better than those 
predicted by Bartz expression. This method has been found 
useful in the present investigation and the details are 
described in Chapter 2. 


i 
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The thermo-chemical response of composite liners 
exposed to some heating rate has been dealt with in various 
ways. Sutton [5] has recently reviewed the initial 
development of ablative thermal protection systems. A few 
experiments conducted under controlled conditions provide 
empirical expressions for the degradation of liners which 
are found useful in design calculations but their range 
of applicability is rather limited. On the other hand, 
detailed investigations involving elaborate computer programs 
are available which attempt to deal with the problem in the 
most general manner possible by considering the non- 
equilibrium chemically-reacting boundary layer flow coupled 
with the chemical-kinetic-controlled heat transfer processes 
in the Charring ablater [6,7]. In view of the complex 
chemistry of the decomposition of composite materials, the 
development of new materials, the lack of accurate knowledge 
of the thermo.-physical properties at the operating conditions 

and large computer time required for 
iterative cycles to achieve coupling, it is useful to 
devise quick and approximate solutions. Goodman [8] employed 
integral method to solve the heat transfer problem in the 
composite heat shield for a given heat flux, and has also 
reviewed other related investigations. 
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The classical problem of heat transfer with phase 
change finds an interesting application in the case, where 
different kinds of techniques have been employed to obtain 
approximate solutions. Newmann [9] obtained an analytic 
solution for a phase change problem with constant surface 
temp Various related problems with specified heat flux 
have been discussed by Boley [10] and more recently by 
Zien [11,12], Here it has been found more convenient to 
use the approximate method of solution, for time-dependent 
surface temperature /by Sharma et al. [13], The details of 
the method are discussed in Chapter 2. 
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CHAPTER 2 

FORMULATION AND METHOD OF SOLUTION 

2.1 PHYSICAL MODEL 

The propellantsin a rocket motor when ignited 
produce hot combustion products which are then accelerated 
through a converging-diverging nozzle to produce thrust. 

After a very short initial transient period, some constant 
stagnation conditions of the gaseous products exist in the 
chamber and the flow through the nozzle remains chocked. 

The flow field in the nozzle is assumed to be steady and is 
analyzed by considering one-dimensional, inviscid 'core* 

(or main) flow along with the turbulent boundary layer 
approximation over the nozzle walls (see Fig. 2). The 
multi-component nature of the flow has been ignored, and 
an approximate method of solution as described in the next 
section has been employed for the turbulent boundary layer 
flow. 

The passage of hot gases over the nozzle wall transmits 
heat into the composite liner. Each composite liner has 
a threshold temperature at which significant decomposition or 
pyrolysis begins to take place. The pre-pyrolysis time, that 
is, the time taken for the liner surface to attain the 
pyrolysis temperature, is usually very low. As the surface 
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temperature of the liner increases further, the temperature 
of the inner layers of the liner rises and therefore pyrolysis 
sets in them, re suiting in the formation of a char layer. The 
gaseous products due to pyrolysis diffuse outwards through 
the char layer, and thereby reducing the heat transfer 
from the hot propellant gases. Assuming that the pyrolysis 
is restricted to a very thin zone, the complex thermo- 
physical processes occurring within the composite liner 
have been described simply as an approximate unsteady heat 
transfer problem with a char front moving into virgin material 
(see Fig, 3) and an approximate method described later has 
been used to solve it. 

2.2 MATHEMATICAL FORMULATION AND ANALYSIS 

* 

The mathematical formulation of the problem consists 
of essentially three main parts, namely, the 'core* flow in 
the nozzle, the compressible turbulent boundary layer flow 
on the nozzle walls, and the heat transfer phenomenon in 
the composite liner. Of course, each is coupled to the 
other through appropriate interface matching conditions. It 
is reasonable to assume a quasi-steady situation in the 
sense that the high speed nozzle flow adjusts itself quickly 
to the changes at the surface of the composite liner in 
which heat transfer process is unsteady. The underlying 
assumptions and the governing equations for each part are 
described below separately. 
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A. ’Core’ Flow : 

The main flow in the nozzle is assumed to be steady, 
one-dimensional and inviscid, and the hot gases are 
Assumed to obey ideal gas equation of state with a suitable 
gas constant as well as to have a constant average specific 
heat. Under these assumptions the governing equations for 
the ’core' flow may be expressed as follows (see, for 
example, Shapiro [14]) : 

Mass conservation ; $ e u e A “ const = m (2,1) 

du dp 

Momentum conservation : 3 e u Q (^~) = - (^~) (2,2) 

Energy conservation i C p T e +r ^ U e “ cons ' t = C p T o (2.3) 


Equation of state : p = 5 

V V t 


(2.4) 


Tho symbols have been defined under the heading 'Nomenclature*. 
These equations can be easily integrated to obtain the 
following results : 


A_ _ r 2_ 

A * ~ M e 1 T+l 


(1+ ^1 M 2 )] ([(Y+1)/2(T - 1)] 


(2.5) 



( 2 . 6 ) 

( 2.7) 


12 


B. Turbulent Boundary Layer Flow i 


For the operating conditions of a typical rocket 
motor, there will exist compressible turbulent boundary 
layer along the nozzle wall. Several approaches have been 
developed over the years for analysing turbulent boundary layer 
flows. In order to avoid the complexity of solving non- 
linear, partial differential form of the turbulent boundary 
layer equations, the so-called ’integral method 1 based on 
the * integral* forms of those equations have been found 
useful especially to predict the wall parameters. The 
integral momentum and energy equations for a compressible 
turbulent boundary layer in axisymmetric flow can be 
easily derived (see, for example, Bartz [l]) and expressed 
respectively as follows i 




*) + i <&> 



( 2 . 8 ) 


@ + > fjt 109 ® 


S V h t,o- h t,w )r 3 


"ad,w 


St (2.9) 


From the equations governing the ’core' flow, eqns* (2.1) to 
(2.4), it can be readily shown that 
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‘ 4 U < 


d(§ u ) 
\r V> e Q* 

;L dx 


] - E 


i-m: 


M g ( 1+ 


Mt )- 


dM 

l(x^) 


ax 


( 2 . 11 ) 


Equations (2.8) and (2.9) can be used to determine the 
Skin friction coefficient (C^.) and the Stanton number (3t) 
provided the velocity and the temperature profiles become 
known. Past experience has shown that the wall parameters 
can be predicted to a reasonable accuracy by using approximate 
profiles such as in the pioneering work due to Bartz [l]„ 

Here, in order to obtain approximate velocity and 
temperature profiles, the more recently acquired knowledge 
in describing the structure of turbulent boundary layer flows 
has been incorporated as suggested by Mastanaiah [4]. The 
earliest models for the phenomenological description of the 
coefficient of momentum transport or the so-called Reynolds 
stress were proposed by Prandtl and Boussinesq in terms 
of mixing length and eddy viscosity parameters. Van Driest 
[lb] extended its applicability over the entire boundary 
layer for flows with negligible pressure gradient and zero 
mass transfer. Cebecci [3] generalized it to flows with 
pressure gradient mass transfer and heat transfer. In 
terms of dimensionless quantities, the generalized eddy 
viscosity expression can be written as 
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where 


D = l-exp(-y + /A + ) and A + = A+( 1-11. 8P + ) 


( 2 . 12 ) 


(2.13) 


In equation (2.13), A* is taken as 26, and in eqn. 
(2.12), K m , the mixing length constant as 0.4 . The tempera- 
ture dependence of the coefficient of viscosity is described 
as T™ ,■ m being equal to 0.6. Here eqn. (2,12) has been 
used to determine the velocity distribution across the 
boundary layer. 

The relationship between temperature and velocity 
profiles in a turbulent boundary layer along supersonic 
nozzle with heat transfer has been investigated by Back and 
Cuffel [16], They have found that the usual Crocco- 
relationship is no longer valid and that the acceleration 
parameter plays a prominent role. Following Mastanaiah [4], 
a quadratic enthalpy-velocity relationship of the form 
given below is assumed : 

^ = A x 9 + A 2 <p ? (2.14) 

where A^ and A^ are determined by using the boundary 
conditions at the wall, namely. 
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At y = o, u - o, T = T w , 



-k( 




(2.15) 


After some algebraic manipulations and also making use of 
the differential form of conservation equations for 
momentum and total enthalpy, it can be shown that 



= 2P r (St/C f ) 'S 


ad,w 




1 - 





+ 



“5 ) 

Q. u h 

-^0 0 o 


d log e (M e ) 
dx 


(2.16) 

The above set of equations, namely equations (2.5), (2.8), 
(2.9), (2.12) and (2.14) are sufficient for the determination 
of skin friction* coefficient and Stanton number for a 
known constant wall temperature. 


C. Gas-Solid Interface Condition : 


The heat flux into the constant- temperature nozzle-wall ■ 
varies along the axis of the nozzle, and at any given 
station it depends upon the magnitude of the wall temperature. 
When the nozzle-wall is lined with a composite material, 
certain thermophysical processes occur at the gas-solid 
interface as described in the last section resulting in 
the coupling of turbulent boundary layer flow with the 
physical phenomena taking place within the composite liner. 
Both the surface temperature and the heat flux change with 
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time soon after the motor is fired. At any instant 
the heat flux from the gaseous phase can be expressed as 

% = h ( T ad -V < 2 - 17 > 


The adiabatic wall temperature can be obtained from the 
relation 



1 + R 

T 0 C ^ -*] 

1 + 


(2.18) 


where, the recovery factor, R is taken equal to . 

At any station on the nozzle wall, T gd is constant while 

* 

q w w iH vary with T w and consequently the magnitude of the 
heat transfer coefficient, h, will also change. The details 
of the coupling Of the heat transfer and the decomposition 
of liner is discussed later. 


When the pyrolysis sets in the composite liner, the 
heat flux at the surface is affected by the gaseous products 
of pyrolysis which diffuse outwards. The hoat blockage due 
to pyrolysis is determined by using a correlation due to 
Fogaroli and Laganelli [17] which takes into account both 
Mach number and wall temperature variations but assumes very 
small pressure gradients. This correlation is as follows i 



[1 



2 , 5 *# 


( 2 . 19 ) 
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where the superscript 'o’ refers to Stanton number without 
pyrolysis-effect and 

b l= ‘S.in/W [<-Sv) w /CSu) e ] (1/St?) 

b 2 = exp[l, 676(0+ 0.l6l)] 

and 

w = I M e + ( V T e )" l/S 

The mass flow rate of the pyrolysis gases, ( $v) w , is 
determined from the kinetic equations describing the 
thermal degradation reactions in the pyrolysis zone. Since 
most polymers degrade in a highly complex manner, semi- 
empirical homogeneous kinetics are normally used to describe 
the degradation by an Arrhenius type of equation. That is, 


»t 


= -A oxp(-E/RT) S„ ( 


< 5,-1 


■)" 


(X.2o) 


D. Heat Transfer in the Liner : 


It is reasonable to assume that the temperature gradient 
in the transverse direction is much greater than that in the 
longitudinal direction. During the pro-pyrolysis period, the 
heat transfer in the composite liner is then governed by the 
unsteady one-dimensional heat conduction equation : 
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, ST, ' 

( 3T-)= 3 V*? ) ' 0<t < V \>° 

subject to the boundary conditions : 

T^XfO) = T° 


( 2 . 21 ) 


( 2 . 22 ) 


T x (d,t) = T° 


(2.23) 


-k 


^ T i 

l S3T 


i- 

X-O 



(2.24) 


A solution of eqns. (2.21), to (2.24) must be 
obtained to determine the pre-pyrolysis time. After the 
onset of pyrolysis, the liner can no longer be treated as 
made up of one substance instead it has a char portion and 
the remaining virgin material such that the thickness of 
char layer increases as the interface moves into the virgin 
material. Assuming constant thermal conductivities (k) 
and thermal diffusivities (3fi») for both the char and the 
virgin material, the heat transfer will be governed by the 
following differential equations. 

t>t , o < x < X (2,25) 


> lh - *. &> 


< 4 £> - *1 4 i> 


3x' 


t > t p . 


X < X < d 


(2.26) 
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subject to the boundary and interface conditions, 
T.(o,t) = T w Ct), t > t p 

Tl(x,tp)= T£(x) 

^(d.t) = T 1 ° 


(2.27) 

(2.28) 

(2.29) 


T c (X,t) = T x (X f t) = T p 


(2.30) 


k 


C3 X 


X 


~ k i $T“ 



(2.31) 


Here, T (t) in eqn. (2.27) is not known explicitly as a 

w 

function of time, t and, in fact, it is responsible for the 
heat transfer coupling between the turbulant boundary layer 
and the composite liner through eqn. (2.17). 


2.3 METHOD OF SOLUTION 


The variables characterising the ’core* flow can be 
readily obtained from eqn. (2.4) to (2.7) provided the area- 
ratio is known at any station along the nozzle. For conical 
section, it can be calculated from the throat radius, the 
length and the convergent/divergent half angles of the 
section, while for contoured section, a polynomial curve, 
fitting the contour is found useful in determining the 
radius at any station. The distribution of Mach number (M e ), 
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pressure (P ), temperature (T ) of ’core’ flow along the 
nozzle is then determined in terms of the chamber stagna~ 
tion conditions P Q and T . 

In order to deal with the turbulent boundary layer 
flow, experimental or assumed values of the quantities & 

©, 'V * and S t are needed at the entrance station 
of the nozzle. Using the values of and af the ith 

-f 4“ 

station of the nozzle, the velocity profile, u versus y , 
can be obtained at the (i+l)th station by solving eqn. 

(2.12) in conjunction with eqns. (2.13) to (2.16). Since 
y + varies more rapidly compared to u + , it Is found conve- 
nient to solve for y for a given step change in u by a 
sixth order Runge-Kutta method (identified as IMSL subroutine 
DVERK). When the value of u approaches 0.9 u Q , smaller 
step size in u + is used to account for large changes in 
y + in this region, and the integration is stopped when u 
attains a value of 0.99 u e . Thus, the velocity and 
temperature profiles are obtained at any station along the 
nozzle. Using the definitions, the displacement, momentum 
and energy thicknesses are then calculated at the (i+l)th 
station by numerical integration. Since the data points 
are arbitrarily spaced, a method due to Gill and Miller [18] 
is employed for integrating function values for arbitrarily 
spaced data points (identified as NAG subroutine D01GAF). 
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Having found the values of 6*, 0 and dp at the 
(i+l)th station, eqn, (2,8) and (2,9) are solved simulta- 
neously by using two-dimensional Newton-Ralphson method 
to get new values of C f and S^, The derivatives, (d©/dsc) 
and d'Vdx) have been approximated by forward finite 
difference expressions so that eqns. (2,8) and (2.9) may 
be expressed as follows : 

& i+l = + £ f i+1 (6 * x ) + f i^ x >] ( 2, 32) 

and 

^i+l = ^i + Cg i+1 (T,x) + gi (%x)3 (2.33) 

Here f and g denote all the terms in the respective equa- 
tions except the derivatives. These values of and 
are now used to calculate new velocity and temperature 
profiles, which, in turn, provide new values of 6*,© and 

t 

UMj ~ 

resulting^new values of and In other words, the 

entire calculation procedure is repeated at the (i+l)th 
station until convergence In the values of and is 
attained with the desired accuracy. 

In the case of a contoured nozzle, it is found 
convenient to change the independent distance co-ordinate 
along the nozzle axis. This can be easily done by making 
use of the following relationship 
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& - [1 + ( 2 - 34 ) 

After having described the method for the determination of 
local heat transfer rate into the nozzle wall for a constant 
wall temperature, the next task is to obtain a solution of 
the heat transfer problem in the composite liner.For this 
purpose, it is assumed that the thermal penetration depth 
is much less than the liner thickness so that the analytic 
solutions for the semi-infinite geometry are employed in 
the analysis. An exact analytic solution of eqns. (2.21) to 
(2.24) for a constant heat flux may then bo obtained as 

given below [9,12]. 

* 

T(x,t) =22 ( ( e -* 2 /4*t _ x u/4iCt) i orfc {x /YTrt)) 

(2.35) 

The above equation provides us with estimation of time 
required to raise the wall temperature from some initial 
value to any higher temperature. It is found from the 
calculations that the pre-pyrolysis time is extremely short, 
much less than one second and therefore negligible compared 
to the run time of the rocket motor. Further, it may also 
be pointed out that since the gaseous flow does not attain 
steady state instantaneously, there exists an overlap of the 
transient period of gaseous flow and the pre-pyrolysis time. 
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It is , therefore, reasonable to neglect it altogether for 
quasi-steady analysis. 

In view of the above analysis, the governing equations 
for char layer formation may be further simplified by 
setting t^'- = o and T|(x) « const = in eqn. (2.28). An 
approximate solution is then obtained as a function of T w (t) 
following the method of Sharma et al. [13] . In this method 
the temperature profile in the char layer Is determined by 
assuming an arbitrary profile which satisfies the governing 
equations, eqn. (2.25)^. and the boundary condition, eqn. 
(2.27)^ at x = o as well as the interface condition, eqn. 
(2.31), at x = X(t) resulting in the expression 

T rT T w(t) i dT ( t) 

T c (x,t> = yt) + [ Jjufr-]* ♦ ^ (-*-> £**<*>]« 

(2.36) 


The temperature gradient at x = o and x = X(t) can then be 
found from eqn. (2.36) as given below 


2L 


rW*), X( t) , 

*■ Wt T J ( 


1 x=o 


2 % 


dt 


(2.37) 


and 



x=oX(t) 


+ *ua (^) 

L x( t) J + 2 Xc ( ^ } 


dt 


(2.38) 


24 


For the virgin material, an approximate estimate for the 
temperature gradient by using the governing equation, 
eqn. (2,26) and the initial condition, eqn, (2,29) is 
found to be 


3 T 

1Tx 


1 


x=X( t) 



(2.39) 


Substituting for the temperature gradients from eqn. (2.38) 
and (2.39) in interface condition, eqn. (2.31) the movement 
of the char front is then given by the equation 


dX£tl 

dt 


dT (t) B[T ( t)-T ] 

A(- H f )X(t) — E. 


JJT 



(2.40) 


where 



c 

2 * c *<5p 




k dV T I> 

?(njt 1 ) 4 ' Q 


(2.41) 


Next, combining oqns. (2.37) and (2.17), the heat-balance 
condition at the gas-char interface can be expressed as 


i dT (t) 
L* w v J 

^dt 


2 3C_ k 2 X 

+ xTtf [ ^% + h]T w (t ) --XTT7 


[hT 


ad 


k; T 

+ 


xltp 


0 


(2.42) 


A solution of the coupled set of ordinary, nonlinear 
differential equations, eqn. (2.40) and eqn. (2.42), will 
give the char layer thickness and the wall temperature as 
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time progresses. For the derivatives, the backward 
difference approximation is employed, that is. 


dJSit 1 \ ■ 

dt | t=t x + a t 


XCt^At) - X(t 1 ) 
“ At 


and 


( 2,43) 


dT w (t) 

dt 


t=t^+ A t 


W At > - W 

__ 


subject to the conditions that 


At t = o, X = o. 


and 



(2.44) 


(2.45) 


In view of the singular nature of eqns. (2.40) and 
(2.42) at t = o, the numerical integrations were started by 
giving small increments to t and T^. Further in the 

derivation of eqn, (2.42), the modification in the heat 

# 

transfer rate, q w , due to pyrolysis-gas blockage as 
expressed by (2.19) is not done because the effect of 
blockage was found negligible. It may be pointed out that 
the dependence of the local heat transfer coefficient on the 
surface temperature is determined as described in Chapter 3, 
and is taken into account while solving the above set of 
coupled equations. 




Fig* 3 Schematic diagram of a charring liner and the co-ordinate 
system used for it* analysis* 
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CHAPTER 3 

RESULTS, DISCUSSION AND CONCLUSIONS 

3,1 RESULT AND DISCUSSIONS 

The computations for predicting heat transfer rates 
have been done for two different situations, namely, flow 
of heated air in a conical nozzle and flow of combustion 
gases in the SLV-3 first stage nozzle. The char layer 
thickness calculations are done for the case of SLV-3 
rocket only. 

Back and Cuff el [ 16 ] measured temperature and velocity 
profiles in a turbulent boundary layer flow of heated air 
in a cooled convergent-divergent nozzle with a half angle of 
10° and with 10:1 contraction ratio. The nozzle throat 
diameter, inlet length and total length are 40.39 mm, 

253 mm and 526 mm respectively. The experiment was conducted 
with P Q = 1.035xl0 6 NnT-2, T q = 835°K and ( T w /T 0 ) = 0.5. 

Under these conditions, the radiative heat transfer is found 
negligible. The Prandtl number is assumed constant and 
taken equal to 0.7. It is claimed that the heat transfer 
measurements are accurate to about + 5 percent and the 
boundary layer thickness to within 10 to 20 percent. The 
comparison of the predicted heat transfer rates and other 
parameters with the above experimental data has been done to 
ascertain the accuracy of the present approximate analysis?. 
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The computer program SARS.FOR is run using the 
experimental data at station 1 as the initial conditions 
with the interval &x = 0.02m along the wall. Computed 
results showing variation of Mach number, displacement, 
momentum and energy thickness© s^kin friction coefficient and 
Stanton number along the axis of the nozzle are plotted in 
Fig. 4 to Fig. 8. The displacement thickness becomes negative 
as seen in Fig, 5 which can be attributed to the cooling 
of the nozzle-wall* The CPU time for one run covering the 
entire nozzle is found around 1.5 minutes using a step size 
of AU + = 2.0 along the Y-direction. 

A comparison of the computed Stanton number due to 
Mastanaiah [4],Bartz [l] (see Eg. (1.1)) and the 
present on» with the experimental data of Back and Cuffel 
[16] is shown in Table 1. 

The- present computed results are 15 percent and about 
6 percent higher than the measured values of Stanton number 
in the convergent and divergent sections of the nozzle res- 
pectively. The results based on Bartz closed-form approxima- 
tion are much higher than the measured values. Improvement 
of present results as compared to those due to Mastanaiah 
may be attributed to the use of more accurate method of solv- 
ing differential equation as well as numerical integration. 
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The comparison of other computed flow parameters 
like 6*, ©, ^ and with those of Back and Cuffel [l6] 
is given in Table 2. 

Unlike the case of Stanton number, the computed skin 
friction coefficient is lower than the experimental values 
in the convergent section. 

After having established the level of accuracy of the 
present method of computing heat transfer rate, the calcula- 
tions were done for the nozzle used in the first stage of 
SLV-3 rocket. The nozzle shape is shown in Fig, 2 having 
throat diameter equal to 330 mm. The stagnation pressure, 
stagnation temperature and the characteristic velocity 
are 3.728xl0 6 Nra” 2 , 3300 °K and 1.5328xl0 3 mS" 1 respectively. 
A cubic polynomial curve-fit is found adequate to describe 
the nozzle contour. Constant average specific heat, 
molecular weight and specific heat ratio of the combustion 
products equal to 1652J Kg” 1 °K” 1 , 25 and 1.25 respectively 
have been used in the calculations. Prandtl number has 
been obtained from the following relationship : 

P r * [4Y/(9*-5)]. • 

. As pointed out in the last chapter, the initial values 
of 6*, are necessary at the first station in order 

to proceed with the solutions of the governing equations 
along the nozzle. No experimental data was available 
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for SLV-3 rocket. However, it was possible to obtain 
solutions with arbitrarily chosen initial values of 
6*, ©,>£ . The trial runs indicated some dependence of 

the final results on the initial guess. Therefore, the 
initial values of 5*, 9, and were chosen such that the 

predicted heat transfer rate at the starting station gave 
values close to the one based on Bartz closed form 
expression. Slight variation of initial values around 
these chosen values had insignificant effect on the heat 
transfer rate. 

The computer program SARS1.FOR is run using the above 
mentioned data for the nozzle-wall temperature equal to the 
decomposition temperature, 823°K of a typical phenolic [19] . 
Computed results showing the variation of Mach number, dis- 
placement, momentum and energy thickness, skin friction- 
coefficient and Stanton number along the axis of the nozzle 
are plotted in Figs. 9 to 13. The displacement thickness 
remains negative throughout the nozzle (see Fig. 10) because 
the chosen value of (T w /T ) is rather low so that the boundary 
layer is highly cooled. The displacement thickness along the 
nozzle ' varies in accordance with the local heat transfer 
rate the variations in the momentum and energy thicknesses 
(see Figs. 11 to 12) are similar and depend upon the local 
heat transfer rate. The skin friction coefficient and 
Stanton number decrease at a rapid rate initially 
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and then increase slightly towards the end of the nozzle 
($ee Fig. 13). 

The local heat transfer coefficient variation along 
the nozzle is shown in Fig. 14 for ^ v /T q = 0.25. The 
maximum occurs slightly upstream of the throat. For the 
sake of comparison the result obtained from Bartz expression 
eqn. (1.1), is also shown in Fig, 14, it may be pointed out 
that Bartz expression has been derived by making the experi- 
mental data at the throat agree with the calculated results. 
The difference in the two predictions for the convergent 
section may be attributed to the lack of well established 
known initial conditions, however, in the divergent portion, 
Bartz expression predicts much higher heat transfer rates 
as compared to the present one. The dependence of heat 
flux rate on wall temperature for heated air data along 
the convergent and divergent sections of the nozzle is 
shown respectively in Figs. 15 and 16. As expected^ the 
higher the wall temperature the lower is the heat flux 
rate and it is nearly linear. 

The char layer thickness calculations were done at 
four different stations along the nozzle which were lined 
with either carbon phenolic or silica phenolic. The physical 
properties and the pyrolysis, kinetic parameters for both 
carbon phenolic and silica phenolic are given in Table 3. 
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A simple computer program CHAR. FOR is run using 
the data given in Table 3 starting with the nozzle wall 
at a certain known temperature. The effect of blocking 
due to the pyrolysis products on the wall heat transfer 
rate [ see eqn. 2„l9)]is found to be less than 1 percent 
and was therefore neglected. 

The results of computations at five different stations 
along the nozzle showing variation of wall temperature and 
char thickness are shown respectively in Figs. 17 and 18. 
Comparison of the calculated char layer thickness with the 
available experimental data is shown in Table 4. The experi- 
mental values quoted within bracket are from reference No. 20, 

close 

where the ^stagnation pressure wtesreh- is^to the one 

in reference 19. But the stagnation temperature is not 
specified, in ref. No. 20. 

The char thickness is found to be more for the carbon 
phenolic as compared to the silica phenolic for nearly the 
same heating rate. The comparison of the calculated thick- 
nesses with the experimental data is better for carbon phenolic 
as compared to the silica phenolic. The variation of char 
depth in silica phenolic towards the aft-end of the nozzle 
is negligible which may be attributed to the nearly constant 
heat flux. 
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3.2 CONCLUSIONS 

Tho general methods of calculating heat transfer rate 
and its effect on the liner materials in a converging- 
diverging nozzle involve extensive use of computer programming/ 
numerical techniques as well as excessive computer time taken 
by the iterative procedure required for matching at the liner 
surface. The present approximate procedures yield reasonably 
accurate results when compared with the available data. The 
method is simple and can be extended to deal with higher 
heating rates involving regression of even the char layer .. and 
for temperature-dependent thermal properties. 
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Data of composite i iners used ln ^ 
computation 

A. Properties of Carbon phenolic [ 19 ] 


* 1270 Kg/m 3 

3 C C = 1067 Kg/m 3 

: . = 1507.2 J/Kg°K 

p, l 

C p>c * 2009.6 J/Kg°K 

« 1,14 J/m-sec°K 

k c = 2,0 J/m~sec°K 

T° =s 330 K 

Tp = 823 K 

: p„inj * 1670 * 5 J / K S° K 

Qp - 77 5x10 3 J/Kg 

Properties of Siller. 

Phenolic [ 19] 

=3 1503.0 Kg/m 3 

3c c = 1327.0 Kg/m 3 

C . « 1256 J/Kg°K 

Pfia 

C » 2009.6 J/Kg°K 

Pf C 

«* 0.5 J/m~sec°K 

k c = 0.7477 J/m-sec°K 

T° as 330 K 

Tp s 823 K 


C p,inJ “ 167D * 5 J/K9 ° K Qp “ 775xl0 3 J/Kg 


C, Thermochemical Data for C-P [24] 


"3S S 5 1 - n 

Equation used ^ = -A exp(-E/RT) S Q (*~ — ~ ) 


S , 1300 Kg/m 3 , S • 770 Kg/m 3 , n = 2 


T < 720K 


, m 

A sec 


T > 720K 
l.lxlO 4 


1 . 1 x 10 


E/R K 


O t 2 

2 # 6xl0 3 
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Table 4 

Comparison of coloured and measured char layer thickness 


Station 

number 

Had. 

mtrs 

Mach 

num- 

ber 

Liner 

mate- 

rial 

Char icxcrfljfc, 
ifeeyear 0 

measured 

Thickness*, cm 
calculated 

1 

0.212 

1.89 

CP** 

1.43 

1.41 

2 

0.230 

2.06 

CP** 

[1.5] 

1.48 

3 

0.26 

2.29 

CP** 

1.2 

1.24 

4 

0.306 

2.57 

SP 

[1.15] 

0.76 

5 

0.318 

2.63 

SP 

1.10 

' 0.76 

6 

0.338 

2.735 

SP 

[1.10] 

\0.759 

7 

0.369 

2.87 

SP 

0.8 

0.69 


CP - carbon phenolic 
SP - silica phenolic 

*Char layer thicknesses were measured at the burn out of the motor 

**At this station* the carbon phenolic is backed by silica 
phenolic whose properties are similar as far as pyrolysis is 
concerned. 
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TFfZ,GT,14.5E-2)fiO TH 1200 
HO TO 1400 

RAD= 19 . 01>>2 + 0.267 95 *fZ»I 4 , 5 e- 2 ) 

CO TO 1 

IF (7* .GT.0,0) WHE=2,!> , _ „ 

RAD*-3l220»§2*Z*Z*Z+1.6tJ0 642*Z*Z~0.35bF>?*Z + l S.bE-7 


ABASTR*CRAD/|AD5TR)fCKAO/RAOSTR) 
S«?.0»(S»MA-I .03/(GAMA+i,0) 


**#* • „ 

IF AREA RAX ID IS LESS. IRAN UNITY THROAT CONDITIONS A 

!SS 53 t!S««B'SJ 03 ii 3 :**S!'a**a;jse»XESS= 3 :s = 3 :q:s = 3 !s:s = 3 :s:s 5 =: 3 ! = *s:as;=-s = =:s! 

Iff ABASTR.LE .1 . 0)QO TO 1 4 

Sl*s -1 ■ 0 ■ 

|*aI|?| 8 ASTR)**S 

• 2 . 0 /(GAMA+ 1 , 0 ) 

fX#<*P#CtMACHE) )**s+5*MACWE*MACHE+COd 
fr*S»MAC«E**Sl + Q*2.0*«ACHE 
:.®|jlil'.*IIAC'RE- ( FX/Ft X ) 
vtftei#««ACHEl-MACHEJ/MACHei 

fficttS*JL8$ e m achf i ) 

IF CABSCFRRDR) ,GT,E) GO TO 11 
TYPE*, Z,At3ASTR,MACHE, RAD 

TF THE AREA RATIO TS GREATER THE PRESCRIBED VALUE ,ST 


Tv f ABASi‘R.nT,6.7) STOP 
0(i l '1 7 07 

»pPuTCAt*I.E n CONICAL NOZZLES ONLY **** 
TtiroAl CONDITIONS ARR GIVEN BELOW **** 



7rR„9 

ALPHJUfl.O 


7 J 7 


l T wITK 


FuPMM’f 'TtlR At *HVF T3 THE MACH NO. A'J' GT W’FT.M 2 ^TA! 


1 ti«H 


90R 

C 


?U NFwTUM PAl.PHSClN MF'l'H JR TS USED TO GET ciFTTgR A PP« 

Bss2:s:£:3:s:s:::srs5;:s:s5sxasssss;B3:ssss:s333=:r3 

continue 

Xff 1)*CFRY? 

xcmasr 

nxsi.uF-*; 

C ALL MYPPGM(XC t nn , s 

Typfe;*,XCUnxcC2),RCfl)»RCCai 

NAA*1 

X»(XC(!)) 

xccn»f l.o-axmetn 

CALL MYPPGW(XC,RCA) 

TyPE*,XCCi)#XC(2) ,PCA(! ) ,RCAf2) 

XCU)*X 


c 

c 




r. 


'.( i ) ) 

WE*,bftX,Of2X,RCA(W,RCAC2) 

jO-DX)*XC;(!J 

?XPg*!?XG 1 1 Yf XC f 25 f RCB 1 1 J ,r:»C2> 

|'2} y 

Sl5gJ«-358iU?^8K?»f«J 

E „ lir tX*OF2T*DP1 XJDP2X 

< 1) * CRC { 2 ) *DF1 Y-.RC C 1) * DE? i ) / DEW 
fCC2)aXC(2) + CRC(l ;*UF2X-.UC.C 21^9^1 X)/r>EM 
WE* „ xc a ) , xc f 2 ) , DF1 Y , DF2 Y 
CALL MYPRGMCXC,RC) 

TYPE*,XCC1) ,XCC2),RCU> ,RC.(2) , DELTA 

TK f lABS(RCf2) J.LF.1 ,OE-6) . AMD . ( ABS ( RC CD ), UF. . 1 .t>K-to) 
1 If f NAA.GT.7) GO XH R1 

KaAbWaA + 1. 

f’u in 7yR 


t-jTuAXTu^;-, C9«-iVE»GFU V* I THAN TRF GiVEN ACCURACY JF SfAi| 


» 



r» JVMOrt /ShAKE/A^TaR, A*. A?,4Xtf,rTP / OFbT,T«,UF,tJrr>rt f ^ 
niM fc Ni>rur ruvci ,Ffii<r’2a‘>t)> i Fu^cUi5^J J .F7£PJ««) 

1 i PkT MF! It) ,C(24J , Wf 2^9) , K(l ) , XCQ) ,RCC*), PCA(? J 
?,PCP12J 

THPE*.riy)f f l),F2X3ff n # «D«rPK,EWH^XK 
PUWLR*^AMA/(^AMA-1 .0) 

ckp k ?« f xr 1 1 j i 
ST*{XC(2)J 
TWJf'l’WtfYirmZEKO 

Tw«bOO.O» AHASTP«f had/0, 16S)**2 

f>^T f T A 5tu 0 £ 

to ftnh’the freestream properties rt ^ 

TEaT7ER0/n ,0+C (GAMA-X. 01/2.0 7 *MACHE*MACH£ ) 

RMACHE*HACRE 

CP**1651 • ftb 

PE*(Te/f.?ERO)**POWEH*PZERO 
tSOR'f ( GAMA*R*TE) J *ma2JH£ 

AIE«*fP*|ZERU 
AXW»CP*fW 
TNOT»TW ; !■ 

TAW«C tiO+Rr*«ACHE*MACHe*(5AMA-l .0)/2.0)*TE 

AI4$«C#*f&R 

TOB Y TlllfciFSS ER Q-/TE 

DkWBDX*? TOBTTE»2.0*fDRBYOZ)/Cf MACHEtMACHE-1 .0)*WAD) 

mtf*CHE.C$.I.O)DL*&f>X*0.0 

DNfiBDXvOLMB&X^MACHE 

0TE8&X*! c TE*( t , 0-GAMA)*MACHE*DMEBDX) /tobw; 
BEliIeAtf-ATW 

-$yg»0j£*{2 . ©*DR8 XDZ*UE)/.CRAD*(MACBE*MACHE-1.0I J 



mwws (TE / TW } * RO W E 
COMM*SORT ( MO 0«T ) * 1 1 .R42E-B 
MUF=C0MM*XE S * £ *0,6 
MUW=C0MM*XW**Q,6 

T ^K*,TZFHn,TW, TF, PE, TEMP, A IE, AIW, A I A W , UMERDX , ME , DLrtB 


’v l n i .. ; 777 

l5X, # TjiEKn,T w ,TE, 

""r.wpriw 

Ti,VftsCP , rtV2*Kn*P , *Pfi*UE 


PE, ATE, AIW,Al A«, 


*75600 r 

*25700 

*25^00 1 

125900 
*25000 
f ?5l 00 
!?b70 ft 
75^00 

?o^ao r 

2560 0 
26600 

26700 r 

?6«00 C 

?69Q0 

2/000 C 

? 7100 C 

?7?00 78 

2/5jO 
27550 

2740 0 

?7500 C 

27600 40 

/7 00 . 

7«00 :.Jfv 

/9y0 
bOOO ■ 

I? ^ g: " ' ~ 

I? a ° J " 

?a > J 

v -7 

.?s»t JO 
170 5 0^ 

? ^ 4 J •"* 


orhw»60K7’(TUW«/KnwWj 

{• TiNjO TtiF VETiuriT/ pRjFiT.ES Tfj rHERUOwOAR X kAYFR 

FUMCH1 ) = i.u 
FUWC?l 1 7*0. U 
FUMC5H 7*0.0 

VI 1 7*0.0 
X»0, 0 
HsO . 5 
Ks 7 

Al*ST»PR*(BiArt-ATW)/Cl7fOTfCFB|2 7 


»2*AF1*AF2+Al 


8«f» 


SilT^R0=fW77(2.0»'i‘W) 


♦ RO*F**i2*V 
)vrCF 8 Y 2 *K 0 t 


p|t 1 R*C £MUF/lowiftnEiuoiB0X7 / ( trrOw **37 
Tfpi*,MUE,5u».R0W£ , UE, 00EBDX. 0TOW. P5TAR . TOW* 

AReMOF.MUW.ROWE.UF.DUFBDX^UTOW') 

■-clirlf*CBdf^uTO«)/Mnw 

ll«|*ll|; ■' 

AJf*f SRlIf 43 v „ 

' "ROW^If»/TgM P 2 * R 0 WW 
Iflf f §»*Fi M'P**0 . 6 

:fiBiOW*fW/TEMP 
!i f yPKU *: *Wm M U W 

TFf (X*0fOW).GE, (0.99*UF) ) OQ TO 100 
T t f tV*HfOrt).GE. (0.9*fJE) 7 ti* 0 „ 1 
V0T.= 1 .OE-4 

”.v = * K . 

t ,r>= i 




GALL OVfcPK f M ,FCN, A,* #XF,»R, Tar., TMR,C,Ww, w,IFr) 
T iPfc;+ # y # n i f j ,o,umppw 


n=x* inn* 

A* fK)*Y(.*! j/CHEFFT 
PHTxM/lfc: 

AiRl*P«Al*p»?|+A?*PHT*PHT 
A A at A T Ft, T + Ai W 

TfcIMp*( Ai.~(f!*«/2.0) )/CP 
-- -y)/(2.0*CPJ 


7X7 

RQ2 


Tw 01 *»fKMP+(U*, - . -_ J[ . ^ 

T*f>E * » AYCK) # U » TEMP# AlfS tR # RHT # TMOT* ROW 
FORMAT ' iris. £7,1# 5X,F6.4,5*,F7 . i 
CUMTTwtJE 

POWafTW/fEMp)*ROWrt 
gUWU«KOw*U 
ROt*EUE*ROWE*UE. 

CUMXPM«ROWU/(ROWF*n£) 
t x *»6 * , go W # « > |D |g W| , OE 


?y» 5 ,, 5 3.u f 0 r c 


IFCCftMfl 


) FUMClCKjsO.O 


rune 2 C K 4>f | if *P|2 * Cost RM 

ruwc?csl*uIo-AisiPj*cgMTRM 


00 


TlfW*ll6W,COW»l»fVtfCUK.5,FUWC2(KJ ,Kl»NC3fK) ,PHT, AlSi'F 
K &K + 1 ' •" ■ f 

^ trifwo'f^e.soo) go to too 

GO TO 404 

97s9ras?ss^ss5:ss3sssssss3sss;s!s:sssjssss 

NUMERICAL INTEGRA IT UN OS TNG NAG ROUTINE RQP ROUNDAR 
INTEGRAL PARAMETERS 

*S5S55Ss;s“52;=;3::ss=:;i:;s;::;r:==:;s:::=;;;::=s:s: 

If AIL * 0 

INTEGRATION FOR DISPLACEMENT THICKNESS 
K ®K->i 

CAT.L DOiGAFCAIrFHNCl ,K, ANS,ER,TFAIT,J 
SUM1=ANS 


integration for MOMENTUM thickness 


TrAj,?. = 0 

D f ,T.u uOiGAPlAK ,F?TNC2,K, ANS, ER, TFA1L) 


II X 


14030 

14100 

!4?0O 


SU M 2 = A«s 


4400 
4S30 
4600 
4700 

mo 
joyo 
)000 
( <jtoo 

15100 
IjlOO 
15600 
3a6yO 
15700 
|15*?00 
15^00 
IblOO 
l6?t> 0 
16300 
1&400 
1t>700 
loOOO 
loOuO 
10 00 
7100 
7? 00 
/ 300 
HOI 
7M 


7b7 

C 

r 

c 


' 


TwTfc.'lKft IT JM Vjp FnFkGY rHirKWtSs 


n jut, * r» 

TYt»F*.fcP 

CaT.l, 001 tafia it # Fn^Ci rKr AM^ER, T{- All, 3 
SUMJ*A^S 


nET.STH«S11«l 

MuM'iHN*SUM? 

FNPGTN*S'»Mr 


TYPE VOgiiSf R.MQMfM.EWRGr*, ROWE, UF. RAD, TAW, Trf 
FURMAX (!?X, , OeLSTK , , 5 xn'*WHK', 5 X, ' FtfRGTK J 
type 757 


Dl?§Oxl<lToYOlPyD|/?pl0lCMACfe»MACHe- 

fUffSS5ft88fWi?8l?2j0?f&$T«/HOMTMK.MACHK*MACHE i*2, 

~ V(M&CHP*WACH§rt .0 0+1.0) 


u,o) n/xroz 


V‘-‘ .^4'-^ 



fO)»gNRGfK*C-ST) 

■»CH* Cffttf-TW 3 /.( XZERO-t W 3 


itrif of ftttfw+r 1 xy c 1 3 +r 1 x y ( 2 ) 0 / xtoz 
mmPM 1 -f 2 m 1 ) + r2 x n 2 > ) / x t a z 


iWssSmMP/DJMTEM 

SrWEWsSfMEW 


FhCt =irFB2M-CFBX2)/CFB2N 
r- j,; <s ,f 55 t Tm V E W - 8 T) / S X W E W 

T »"i0+.r r ni? f CFB2 f USX,STNEW # F1 XY( 1 ) ,FtXY(23 ,F2XY(l » ,*"2 
Y-xPb*,Sl ,L»THKTA,f>pSi ,L/BYU|)X,DIMTEM,STTFMP 
i) = ( CFt*?<»)-rFRY2 


\|F«)-ST 

’ r ')OE*,nmE^A,oPST ,fktf, erst, si, dbyitdx 


t* II 


PhTyPrt 

r„n 


rsassssxsssssssissstsrsssxsssBssssssssasssssssasliss 

; 3 ::s:srsE;:sE 3 : 3 Bi 3 :ssss 3 :x 3 : 3 a»s 8 ::as 333 ssmsS 3 s 

SUBROUTINE FCNCN.X.ytXPRlMS) 

COMMON / 5 HftRK/ASTAR f fli # A 2 *AIW f cp # uFLT,TW#UK»ffrnrt- 
DIMENSION V ( N J , IfPRT MF( N) 

AAAX«X*UTOW 

PHT»AAAX/UK 

AIRTAR*A1 ♦PHI+A2*PHI*PHI 
A1»AT5TAR*PELI4-AI W 
TEMpBCAI^AAAXiAAAX/a.Ol/CP 

B*rfEMP/TiO**l,2 


i 


1 


SUBFNsSOR ftg'M.U* HO0RU W*0. 1 6** f N ) *Y( NJ *OMPFN *DMpF.v ) 

■tpf»tinttijll!g 6 ;**'U 0 *T{B) +sa 0 FN) 





*?* ?? , 5fffiR#?S»3555RfJ5S^S^S4Sj5S35JS55'i5SS5JRSSJ5Sri!SS55 
552 5 ^ 5555535555555555555555 ; 

S^£iS 38 SS5S ASSESS 5525 S5 5 


5 55 S 55 S 55 S’ *■* S 55 S 55 
38S*583*55®53SBJS 
38 855 888^S85SS 


i 


»sss*s = sss*a*asaS8Ss*saaa!ss ii*t*i*«**«tst*f*********| 

chargor I 

a***tassi:s*ssi*8S«i***t«a*s**BS«i#ii***»«i«®*»«**** #I j 

T H T6 PROGRAM CAL.rut.ATKS Trie CHAR LATER THICKNESS j 

as wrut. as the surface temp nr acharring ablator 

*aa»*s«s|ca»*a»»a«a**«*®»*a«*a***«****«*"*f*** , M»*u*** 

Tufr subroutine ARR calculates the rate at mrtc« ok; 

subroutine arrr^calchla tes the slocking effect | 

oimenston att( io) . oRosirrc 30) I 

C,lVE n TH? n iIpttT OA£A FROM MEAT TRANSFER ANALYSIS ; 


, IT NF«3jot 0 | 
R23.0?V£R<«t 


CPVIRa 
CHARM 

0bps7 I ^VVV* V. * JWV« VI 

«. <m m mm m f* mm W*rt» • m -» W «•<■- ">"» •* 

T I M EP 0 * 0 f 0T I HE B 0 * 0 0 0 0 1 

0 ZERO | * AH* t T2ER0 *• f I N E ) 


§m& 3 J H«Q*£ P , 

M8«Ai»eii4iwatR:qjf¥lR*CRCHp) 
mwm-2WTm-i$ §g»vi **££?! R) 


?ISIf|RilStTf-II»F)/(0ZER6a*1.12S*S0RT()U,PHBl)) 


f;' 1 -’. !.;;V ‘j »4 C 
J ■-> Cr.'<-. r i 

&MM 



m«*£TP-TINF3/£0ZERh*1 .1 28*SQRT( AT.PHAt )) 

aim*. ■ 


Wititmt L TIHE2 , TIME 

M m'i** of I ME«0. 05 ; CHARDoO . 0 

AAaCrtA«|C?£2.0*0LP*RO«VlR*ALPRA1 ) 

RtiaCHARK/COLP*ROWVTR) , 

r'C=O?IPK*(TP-TINF))/fQLP*R0NVIR*SURTf3.1416*ALPHA?) > 


U „ T> | 


GT VR T rlE DATA FRQH THE NEAT TRANSFER PROGRAM 

_ \ r r* * * t A -» /*» r\ fi" 


A ’iJ *>170 = 2465. 0» Ari.R2lBl438.0S 


7 .< 7 


n f»uTsfl3ft.n; >\CHaB 0 = 1 .OE-B 
' r .l, , V,= j'TMP+nTTMB' 

TffATNOv .00.810,0)30 TO 91 



WQO 

i|«jO 

'50 
6100 
5?0O 
5)00 
54^0 
'5500 
5600 
5^00 
6000 
6000 
ion 

'6500 
! ofiOO 
60 GO 
117000 
07130 
f^on 
07)30 
07400 
07530 
1/730 
)/o on 
1/000 
'6000 
H 1 30 
o530 

3630 

a^oo 

«030 
oOOQ - 
30||, ; 

|qr)o 

; " i O 


3 > -J n 
3 0 7 " 
J s j ’ 

U 6 ,! 0 

7 ) O 

V V* 

t 0 , 1 o 




1 20 
1 2^4 


* t*uiTsT nh r*20 .0 } acharu*Charo+i .of-** 

Tt- f rTrtE.cT,o.o*)r>rTHFaim*F*to.Q 

ifitime gM.i iimuFsO.os 

TFfiTMF.CT.l.O)0TlMe«0,5 
TFf TTMF.nr. 8.0) DTI«fe;*l ■ 0 

pAKTlOaRH^I^*ArHftHr>+?rtRRK»Tp + TNOT*ACHARU*ArHARD/f 2 
1 *CHARK/rvrTMF 

PARTI 0«<AHS23 + {TK3T-82 3.0) *UH?17Q-AH823)/UB00, 0-82 
tArHRRO*Tl3nr/(2.0lAbPKBf*0T!»<e>+CH»R»f*TP/ACHAR0 u ^ 
PARI) 1 »CH ARK/AT HARO+ AC H AR0/( 2,0* ALPHA 1 *DTIM£)+(ARh 2^ 
1 873. 0)*fAH2170-AH823)/< I 800. 0-823.O)) 
ATWOfaPARriO/PARTli 

AAA»f ) ,0/OTIME+AA* C ATNOT-THOT) /DTI Mg) 
PBP*CC£/S0RttTlMF)-CHAR0/OTXPe) 

CCC»fcM8#tAf»&f-rP) 


'g At f t tSif 4 *f * AS A* CCC. 

ACHAR0bC?#rI+|PRT(PAR|)2O)/(2.*AAA) 

IfHIIItlfliRO-OLOt )/IcHAR 0 ).LT. I. OB-6) GO 10 R88 
OLD1|?ACH1RO?OLD2 = AT|SOT 
. 0n xff"' '99-' 

PAM.# f A€ H ARP -CH A RO ) / DTI M£ 
t»|if*AfNt)T?CHARO*ACHARD 

j«>OT*tAHf23+CTMOT-823. ) * ( AH2) 7 0-AH823 ) / ( 1 800 . -R23 . ) )* 
l|PE. CHARD, Tfiof 

GO TO *07 
rOMTTN0£ r . 

FURMPTfS*. 'BLOCKING EFFECT IS CALCULATED BELOrf') 

C A T, L» ARRf CHARD, ALPHA2, RAP, R0VW) 

RQW£HE”92i , 0. ;TW»2170. 1 A.«ACHjri 

CALL ARRR * ROVW, ROWEfJE » TW, AH » APACHE, AH,1 ) 

c/jnp 

r-'.-iH 

:: = 3;:::ss;:;s::~s:s33SSZ5sss3i3ias33SS3ssBS3Ssi>a?8szi 

*'A c -b ET.ijW 1 Jr TO EHE ft 3LATIPG SHREACE 


a * rs*s» as =**» s* a* **c* s***«*8 a****************** **** 3 


SunKnu^lNe ARKfCHARD.ALPHAT.RAMfROVirf) I 

n i M t:N 6 TuM n* Y( 5 u) , DRnaYrf * 0 ) j 

TpT!sl'o7^oIjPa690.0?Tfc,**P«rPii + 20.0| W«1 f T*>ASr.*2| 7 0 *0 t fj 
ChAR* *1 * * | 

a dfl* RAM/ ALPHA? ; 

I^A*un( rr?M^«Ti«F)/rrPAftb*rTNF))/»oFL ! 


*¥YCN)*CHAP0-AV 


TFf TEMP ,CT*720.0) GO TO. I 
TFCTEMP . LC* 720.0) GO 10 2 

rn | 


>„(MIN0EX*2 

VfixnG ] i-mniim variation or density with temp from v 


RO«*fROwR*ROZeRO> * ( TEMP -.TP) / ( TPL-TP) TROZERO 

#St01YT<N)^ASSC-A*SXF(-Ei:BYR/irMP)*ROZF:Rn*( (ROW-ROWR)/r; 

lint )--' - 


■ill '™ IS TAKEN BECAUSE THE DENSITY TS DECPE.ASIN 

fc(TEMP.LE.TP)GO TO 4 

‘ “‘"•5 


?IM r 'DofGAF(AYY,DRn8Y'f , NEW, ANS,ER, IFA ID 


continue 


Tb CHAR REMOVAL IS TAKING PLACE USE THE EQUATION RfL 

^ub*R=07ERQ/QLABL* (ANEWC*C3. 0-2. 0/DEL)/ U.OTANE.WC)) 

PuVw=AHS+REMCHR 


O 


Pfc;T U PiM ! 

j 


s mn 


?STS“S-?ffif“S 3 S 5 B 3 ?I 5 a?fSS"?I 3 SSSS"B 3 r?K-SCSsjEW: 

r"rgrT one r?) thr pfmn,KSIs cases 


SUBRHUTIHe A&RR (ftOVH r ROUfRUR# Trt * AH, AMACHE* AHl 5 
CPT*jIi^8.5FCPAtR«1651.«6>STHhr*AH/CROWRUR»CPA|Rl 

IIfei!^* MA=,,RfwCHe, 

W$W9 g MM % Z A S- Sm+ZJLW n 


It? 




OG5ft5sST|SR*.aO»eUE» C rAtf-Tt&RO) 




■' 'r’ft 

^i M p^xty «saaiegsBS'asgggB>-g-- B - g »* wa ..agi^»y ***■’«■* T’-ISZZ^J 

•^ a ^^«f^#.t#«^.W*|aa***s*** a *»*a**** 6 *****5 s » s ** !Sa * :? *** i, ‘*‘** , 1 


ssa^assssscs^srssssassfsssssBSBsssssfersa 


